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I.  INTRODUCTION 


The  requirement  for  rational  analytic  procedures  to  predict  damage  to  lightly 
armored  vehicles  and  aircraft  caused  by  high  explosive  warhead  detonations  in  close 
proximity  to  such  targets  provided  the  motivation  for  the  investigation  to  be  reported. 
It  would  obviously  simplify  the  development  of  such  procedures  if  the  loading  produced 
by  the  warhead  could  be  predicted  independently  from  the  structural  response  analysis. 
While  this  uncoupling  cannot  be  justified  in  general,  the  extremely  short  duration  of  the 
positive  phase  of  blast  pulses  from  small  explosive  charges  in  comparison  to  the  time  for 
r  opreciable  response  of  moderately  hard  targets  suggests  that  this  simplification  may  be 
suitable  for  the  contemplated  application.  However,  it  must  be  recognized  that  even 
when  the  blast  pulse  approaches  a  pure  impulse,  the  relative  delivery  times  of  such 
pulses  at  different  portions  of  the  structure  must  be  adequately  represented. 

While  it  is  recognized  that  fragments  from  detonating  warheads  are  effective  in  pro¬ 
ducing  target  damage  and  interact  synergistically  with  airblast  effects  on  targets,  it  was 
preferred  to  start  this  investigation  by  attempting  to  define  the  load  distribution  pro¬ 
duced  by  bare  charges  since  it  was  this  augmenting  damage  mechanism  which  had  not 
been  taken  into  account  in  earlier  vulnerability  models.  Also,  there  is  interest  in  the 
effects  of  lightly  cased  munitions  for  which  the  contribution  of  fragments  to  the  total 
loading  is  insignificant.  Consequently,  attention  was  directed  to  the  idealized  problem 
of  defining  the  blast  loading  on  a  large  flat  surface  produced  by  center-initiated  bare 
spherical  charges  at  various  stand-offs.  It  had  been  expected  that  it  would  be  possible 
to  produce  an  empirically-based  computer  program  which  would  provide  the  blast 
pressure-time  history  for  given  input  parameters.  In  fact,  Mr.  K.  O.  Opalka  has 
developed  such  a  program,  but  its  release  has  been  deferred  owing  to  doubts  regarding 
its  data  base  for  the  region  around  the  transition  from  regular  to  Mach  reflection  and 
for  small  scaled  heights-of-burst.  The  existing  data  gaps  and  contradictions  will  be  dis¬ 
cussed  in  the  next  section. 

It  will  be  seen  that  prediction  of  normally  reflected  pressure  and  impulse  for  small 
Z  is  clouded  in  considerable  uncertainty,  i.e.,  there  are  large  differences  between  scaled 
experimental  data  and  a  corresponding  hydrocode  prediction.  This  prompted  an  effort 
to  obtain  blast  data  for  the  cited  small  Z  range.  However,  the  present  status  of  pressure 
transducers  appeared  to  preclude  measurements  in  this  very  high  pressure  region.  On 


*  Hopkinson1  scaled  distances  Z  less  than  0.2  m/kg1/3  (0.5  ft/lb1/3). 
1  B.  Hopkinson,  British  Ordnance  Board  Minutes  13565,1915. 
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the  other  hand,  it  seemed  entirely  feasible  to  employ  a  previously  used  impulse  plug 
technique^  to  obtain  data  for  reflected  impulse.  Such  experiments  were  performed  and 
the  resulting  information  and  comparisons  with  other  sources  are  given  in  the  sequel. 


H.  REVIEW  OF  EXISTING  INFORMATION 

Q 

An  examination  of  the  widely  referenced  compilation  of  blast  data  by  Goodman 
reveals  that  none  of  these  data  involved  actual  pressure  or  impulse  measurements  for 
Z  <  0.2  m/kg1'  (0.5  ft/ib* '  ).  Certain  information  on  normally  reflected  impulse  for 
smaller  Z  is  listed  but  reference  to  the  source4  reveals  that  this  information  was  based 
on  Sachs’  scaling0  of  data  obtained  in  a  reduced  pressure  environment  and  no  shots 
were  fired  closer  than  the  Z  cited  above.  Subsequent  data  from  tests. at  BRL  appear  to 
have  the  same  limitation. 

Some  closer  range  pressure  and  impulse  data  have  been  reported0  by  Southwest 
Research  Institute  personnel.  These  data  provided  information  on  reflected  peak  pres¬ 
sure  and  reflected  impulse  at  various  radial  distances  from  “ground  zero”  along  a  flat 
surface.  The  closest  reported  scaled  distance  was  0.12  m/kg1/3  (O.Sft/lb1/3). 

Another  source  of  information  for  reflected  pressures  is  the  report  by  Kingery  and 
Pannill/  which  is  based  on  the  theoretical  treatment  of  regular  reflection  by  J.  von 
Neumann.  This  report  provides  tables  of  incident  and  reflected  overpressure,  dynamic 


O.  T.  Johnson,  J.  D.  Patterson  II,  and  W.  C.  Olson,  “A  Simple  Mechanical  Method 
for  Measuring  the  Reflected  Impulse  of  Air  Blast  Waves,”  Ballistic  Research 
Laboratories  Memorandum  Report  No.  1088,  July  1957. 

3 

H.  J.  Goodman,  “Compiled  Free-Air  Blast  Data  on  Bare  Spherical  Pentolite,”  Ballistic 
Research  Laboratories  Report  No.  1092,  February  1960. 

4  W.  C.  Olson,  J.  D.  Patterson  II,  and  J.  S.  Williams,  “The  Effect  of  Atmospheric 
Pressure  on  the  Reflected  Impulse  from  Air  Blast  Waves,”  Ballistic  Research 
Laboratories  Memorandum  Report  No.  1241,  January  1960. 

^  R.  G.  Sachs,  “The  Dependence  of  Blast  on  Ambient  Pressure  and  Temperature,” 
Ballistic  Research  Laboratory  Report  No.  466,  May  1944. 

6 

J.  J.  Kulesz,  E.  D.  Esparza,  and  A.  B.  Wenzel,  “Blast  Measurements  at  Close  Standoff 
Distances  for  Various  Explosive  Geometries,”  Minutes  of  the  Eighteenth  Explosives 
Safety  Seminar,  Vol.  I,  pp.  405-445,  September  1978. 

v 

C.  N.  Kingery  and  B.  F.  Pannill,  “Parametric  Analysis  of  the  Regular  Reflection  of 
Air  Blast,"  Ballistic  Research  Laboratories  Report  No.  1249,  June  1964. 

g 

R.  Courant  and  K.  O.  Friedrichs,  “Supersonic  Flow  and  Shock  Waves,”  Interscience 
Publishers,  New  York,  1948,  pp.  327-331. 
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pressure  and  angles  of  the  incident  and  reflected  shocks  as  a  function  of  horizontal  dis¬ 
tance  for  selected  scaled  heights-of-burst.  These  values  were  computed  using  the  von 
Neumann  model  and  data  for  7  (ratio  of  specific  heats)  for  air  as  a  function  of  overpres¬ 
sure.  The  model  does  not  take  account  of  explosion  products  so  the  results  may  not  be 
representative  of  an  actual  explosion.  - 

The  writer  felt  it  would  be  useful  to  make  predictions  of  reflected  pressure  and 
impulse  for  a  small  Z  appropriate  for  a  warhead  detonation,  employing  all  available 
sources  and  extrapolations  from  these.  The  case  of  blast  from  a  0.907  kg  (2  pound) 
50/50  Pentolite  charge  detonated  with  its  center  at  height  H  of  63.5  mm  (2.5inches) 
above  an  infinite  rigid  plane  was  chosen  for  comparisons. 

At  the  writer’s  request  Mr.  R.  E.  Lottero  employed  the  HULL®  hydrocode  to  obtain 
a  numerical  solution  for  this  case.  The  results  of  his  calculations  are  compared  with 
predictions  based  on  the  previously  cited  sources  (or  extrapolations  therefrom)  in  the  fol¬ 
lowing  figures.  Figure  1  shows  the  peak  reflected  pressure  on  the  plane  surface  as  a 
function  of  the  radial  distance  along  the  plane  from  “ground  zero.”  The  corresponding 
reflected  impulse  distribution  is  represented  in  Figure  2.  The  most  significant  results  of 
these  comparisons  is  that  the  peak  normally  reflected  pressure  predicted  by  the  hydro¬ 
code  is  greater  than  any  of  the  empirically  derived  values  while  for  the  normally 
reflected  impulse  the  converse  is  true,  the  differences  being  as  much  as  a  factor  of  two. 
The  figures  also  show  that  these  discrepancies  exist  not  only  directly  under  the  charge 
but  also  over  much  of  the  range  of  R/H  where  the  pressures  and  impulses  are  large. 
Consequently,  until  these  discrepancies  can  be  resolved,  structural  response  predictions 
based  on  such  distributions  are  also  subject  to  serious  doubts. 

Before  arriving  at  conclusions  regarding  the  significance  of  the  variance  between  the 
curves  in  Figures  1  and  2,  it  is  appropriate  to  take  a  critical  look  at  the  basis  for  each  of 
these.  Also,  there  are  three  fluid  dynamic  regimes  to  consider: 

(1)  The  region  of  regular  shock  reflection  which  applies  until  the  angle  at 
which  the  incident  shock  strikes  the  plane  reaches  a  critical  value,  for  this 
case  corresponding  to  R/H  ^  0.9. 

(2)  The  region  of  irregular  Mach  reflection  which  features  a  two  triple  point 
shock  structure.  The  upper  limit  on  R/H  for  this  region  is  not  known  for 
this  case  but  is  reached  when  the  two  triple  points  coalesce. 

(3)  Beyond  the  irregular  Mach  region  there  is  the  regular  Mach  reflection  zone 
which  applies  for  all  greater  values  of  R/H. 


9  D.  A.  Matuska  and  R.  E.  Durrett,  “The  HULL  Code,  A  Finite  Difference  Solution  to 
the  Equations  of  Continuum  Mechanics,”  AFATL-TR-78-125,  November  1978. 


The  curve  in  Figure  1  labeled  Kingery  and  Pannill  was  obtained  by  interpolation  of 
their  tables,  which  are  based  on  von  Neumann’s  theory  for  the  regular  reflection  regime. 
These  tables  apply  to  shocks  in  air,  with  the  ratio  of  specific  heats  7  treated  as  a  func¬ 
tion  of  incident  overpressure.  Since,  for  small  Z,  the  explosion  products  reach  the  plane 
surface  it  is  believed  that  the  reflected  overpressures  are  actually  somewhat  greater. 
The  von  Neumann  model  does  not  consider  the  tail  of  the  blast  wave  so  no  information 
on  reflected  impulse  can  be  derived  thereby,. 

The  curves  representing  extrapolations  of  Southwest  Research  Institute  data  were 
provided  by  Mr.  C.  King^ry.  While  these  extrapolations  were  carefully  performed,  the 
writer  has  determined  through  examination  of  the  scatter  in  the  original  data  and  tak¬ 
ing  account  of  the  fact  that  the  “nearest”  experimental  data  are  for 


REFLECTED  IMPULSE 


Z  =  0.12m/kg1/3  (O.Sft/lb1/3)  that  a  broad  band  centered  on  these  curves  would  be 
more  appropriate. 

The  single  data  point  oh  Figure. 2  at  R/H  =  0  was  obtained  by  interpolation  using 
Table  B  of  a  report  by  Soroka.10  However,  the  data  source  for  this  portion  of  the  table 


*  B.  Soroka,  “Air  Blast  Tables  for  Spherical  50/50  Pentolite  Charges  at  Side-On  and 
Normal  Incidence,”  Ballistic  Research  Laboratory  Memorandum  Report  ARBRL- 
MR  02975,  December  1979. 


is  the  previously  cited  report  by  Olson,  Patterson,  and  Williams  concerned  with 
impulse  plug  experiments  performed  in  an  evacuated  chamber.  The  authors  of  this 
report  were  rightly  concerned  about  a  possible  misinterpretation  of  their  data,  for  if  sea 
level  atmospheric  pressure  is  substituted  in  Sachs’  scaling  parameter  one  might  infer 
that  their  data  included  a  firing  where  the  height-of-burst  was  less  than  the  radius  of 
the  explosive  charge! 

For  large  values  of  R/H  (i.e.,  well  into  the  regular  Mach  reflection  regime)  the 
HULL  code  predictions  and  the  extrapolated  Southwest  Research  Institute  data  tend  to 
coalesce.  This  should  occur  because  at  these  ranges  accurate  experimental  measure¬ 
ments  of  pressures  can  be  made  with  piezoelectric  transducers  and  because  hydrocode 
solutions  become  insensitive  to  the  modeling  of  the  initiation  process.  More 
significantly,  these  predictions  appear  to  asymptotically  approach  values  associated  with 
the  limiting  situation  depicted  in  Figure  3.  In  Figure  3(a)  the  shock  structure  above  an 
ideal  reflecting  plane  for  the  regular  Mach  reflection  region  is  shown  for  a  spherical 
charge  of  weight  W  at  a  low  height-of-burst  H.  Obviously  the  same  shock  structure 
would  occur  if  the  reflecting  plane  were  removed  and  a  second  charge  of  weight  W  and 
height-of-burst  -H  was  detonated  simultaneously.  At  sufficiently  large  values  of  R  the 
blast  parameters  for  the  case  (Figure  3(a))  shown  become  the  same  as  those  for  the  case 
shown  in  Figure  3(b);  i.e.,  for  a  free  field  burst  of  charge  weight  2W  at  the  point  0. 
This  suggests  that,  for  large  R/H,  the  blast  loading  on  the  panel  surface  can  be  con¬ 
structed  using  tabulated  side-on  blast  parameters  for  a  doubled  charge  weight. 


INCIDENT  SHOCK 
REFLECTED  SHOCK 


Figure  3.  Equivalent  Air  Blasts  at  Large  R 


Mr.  Opalka’s  code  attempted  to  combine  various  empirical  models  to  provide  load 
predictions  for  the  three  shock  reflection  regimes  previously  cited.  For  the  regular 
reflection  region  his  code  employed  the  von  Newmann  model  for  peak  reflected  pressures 
and,  using  the  same  spatial  distribution,  determined  impulses  to  match  tabulated  data 
for  reflected  normally  impulse  (which  assumed  the  Sachs’  scaling  law).  For  the  Mach 
region  the  code  employed  an  adaptation  of  an  empirical  model  due  to  Moore,  which 
was  not  derived  for  blast  from  spherical  charges  and  which  seems  to  overpredict  the 
loading  for  large  R/H. 

In  summary,  the  writer  found  no  experimental  data  which  were  obtained  for  the 
small  scaled  distance  range  of  interest.  Extrapolation  from  data  for  larger  scaled  dis¬ 
tances  or  acceptance  of  Sachs’  scaling  lead  to  large  differences  from  hydrocode  predic¬ 
tions  in  both  peak  reflected  pressure  and  reflected  impulse.  The  modest  test  program 
presented  in  the  next  section  was  undertaken  for  the  purpose  of  resolving  at  least  some 
of  these  discrepancies. 


ffl.  THE  EXPERIMENTAL  PROGRAM 


A.  Test  Facility 

Since  pressure  transducers  are  not  available  which  can  survive  and  function  in  the 
very  high  pressure  region  directly  under  the  charge  (cf,  Figure  1)  it  was  decided  to 
employ  the  impulse  ^lug  technique  in  a  similar  manner  to  that  reported  by  Johnson, 
Patterson  and  Olson/  To  avoid  problems  with  plugs  binding  in  holes  the  method  was 
used  only  to  measure  normally  reflected  impulse.  The  analysis  employed  to  interpret 
the  test  data  will  be  presented  later  in  this  report.  The  facility  which  had  been  built  for 
the  previous  tests  was  still  available  and  was  modified  for  use  in  this  series  (see  Figure 
4).  This  facility  was  basically  a  cubicle  with  a  roof  of  25.4  mm  (1  inch)  thick  armor 
plate  containing  a  cylindrical  hole  through  which  the  plug  was  projected  by  the  explo¬ 
sive  blast.  Some  preliminary  firings  were  performed  for  Z  =  0.2  m/kg*'3  (0.5  ft/lb*'3) 
in  the  original  facility  but  it  was  realized  that  the  roof  would  be  inadequate  for  smaller 
scaled  distances  so  most  of  the  roof  was  replaced  by  a  slab  of  76.2  mm  (3  inch)  rolled 
homogeneous  steel  armor  plate.  This  slab  was  of  sufficient  size  that  the  effects  of  shock 
diffraction  at  the  edge  would  not  have  time  to  propagate  back  to  the  plug  during  the 
positive  phase  of  the  air  blast. 

Both  Figures  4  and  5  show  the  metal  frame  used  to  support  the  explosive  charge. 
In  figure  5  an  explosive  charge  has  been  taped  to  a  cardboard  tube  which  is  affixed  to 
the  vertical  adjustment  rod  of  the  support  frame.  Since  the  standoff  becomes  very  criti¬ 
cal  for  small  Z  a  special  adjustable  feeler  gage  was  developed  to  control  the  gap  between 
the  plug  and  the  surface  of  the  spherical  charge. 


**G.  R.  Moore,  “Calculations  of  the  Reflected  Overpressure  and  Transient  Loading  on  a 
Deck  from  Elevated  Gun  Fire,”  Naval  Weapons  Laboratory  Technical  Report  TR- 
2847,  November  1972. 
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Figure  6.  Schematic  Section  Through  Test  Facility 


Figure  6  is  a  schematic  drawing  of  the  interior  of  the  cubicle.  A  nominal  12.7  mm 
(0.500  inch)  diameter  hole  was  bored  through  the  roof,  then  a  15  °  conical  reamer  was 
used  as  indicated  to  minimize  interference  with  the  plug  during  projection.  The  plugs 
were  recovered  in  a  box  on  the  floor  of  the  cubicle  which  was  filled  with  celotex  and  ply¬ 
wood. 

B.  Instrumentation 

It  was  preferred  to  determine  the  plug  velocity  by  purely  optical  means.  Accord¬ 
ingly,  a  vertical  scale  was  mounted  on  a  board  closely  paralleling,  the  expected  path  of 
the  plug.  Two  capacitance  discharge  (spark  gap)  flash  units  were  positioned  as 


indicated  in  Figure  6  and  a  Polaroid  camera  with  a  Wollensack  162  mm/F4.5  lens  was 
mounted  so  as  to  view  the  plug  passing  in  front  of  the  scale.  During  a  test  the  cubicle  is 
darkened  and  an  automatic  sequencer  opens  the  camera  shutter  as  the  firing  begins, 
causes  each  flash  unit  to  function  as  the  plug  is  at  about  its  level  and  then  closes  the 
camera  shutter.  The  precise  time  of  each  flash  is  recorded  on  a  Nicolet  digital  oscillo¬ 
scope.  Since  the  flash  durations  are  sufficiently  short  to  “stop”  the  plug  motion  the 
corresponding  positions  of  the  plug  can  be  read  from  the  Polaroid  print  which  include 
the  vertical  scale. 

C.  The  Impulse  Plugs 

The  following  cons,  erations  were  taken  into  account  in  the  design  of  the  plug 
specimens: 

1.  The  pressure  distribution  on  the  top  of  the  plug  should  be  as  uniform  as  possible; 
this  dictated  reducing  the  diameter  of  the  plug  to  about  12.7  mm  (0.500  inch)  diame 
ter  due  to  the  close  proximity  of  the  spherical  charge. 

2.  The  clearance  between  plug  and  hole  should  be  minimized  to  avoid  diffractive 
unloading  of  the  top  surface  of  the  plug  and  blow-by  of  explosion  products,  yet  be 
sufficient  that  Poisson-type  deformation  of  plate  and  plug  under  stress  should  not 
produce  frictional  resistance  to  motion  of  the  plug.  As  a  compromise  the  following 
machining  dimensions  were  selected: 

Plug  diameter  12.45  +  ♦  mm  (0.490  +  *  in) 

— • UZD  " • UUl 

Hole  diameter  12.70  +_  *  qoo  nun  (0.500  * '  qoo  *n) 

For  elastic  deformation  including  dynamic  overshoot  these  tolerances  should  provide 
0.127  mm  (0.005  in)  clearance  for  a  concentrically  positioned  plug. 

3.  The  plug  should  remain  intact  and  not  experience  sign ifr ant  plastic  deformation 
when  subjected  to  longitudinal  stress  waves  resulting  from  blast  loading.  Initially, 
rolled  homogeneous  armor  was  selected  as  the  plug  material  on  the  basis  of  its 
toughness.  These  plugs  performed  satisfactorily  for  the  larger  scaled  distances  in  this 
series  but  as  Z  was  reduced  it  was  observed  that  appreciable  radial  plastic  deforma¬ 
tion  occurred  near  the  loaded  end.  Following  this,  plugs  were  machined  from  Bear¬ 
cat  tool  steel  which  has  a  higher  yield  strength.  No  problems  with  either  plastic 
deformation  or  fracture  were  experienced  with  the  latter  specimens. 

4.  The  plug  should  acquire  a  velocity  within  a  range  of  values  which  can  be  determined 
with  sufficient  accuracy  by  the  available  apparatus.  The  plug  velocity  depends  upon 
the  blast  impulse  and  its  mass.  The  latter  can  be  controlled  by  varying  the  length  of 
the  plug;  it  was  found  that  use  of  25.4  mm  (1  inch)  long  plugs  produced  satisfactory 
velocities  for  the  entire  range  of  scaled  distances  tested. 
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The  analysis  required  for  interpretation  of  the  test  data  is  considered  in  detail  in 

Appendix  A.  By  definition  the  impulse  delivered  to  the  plug  (actually,  impulse  per  unit 

area)  is: 

I  =  /0Tp^dt 

This  impulse  can  be  evaluated  for  individual  test  data  sets  by  use  of 

(2) 

provided  that  the  following  assumptions  are  satisfied  for  the  particular  test; 

1.  The  displacement  of  the  plug  during  the  blast  pulse  must  be  very  small.  This 
assumption  can  be  evaluated  for  specific  tests  by  using  the  formulation  given  in 
Appendix  A  with  plausible  values  of  the  pulse  duration. 

2.  Friction  between  plug  and  plate  are  negligible.  This  is  necessary  to  obtain  satisfac¬ 
tory  values  for  the  impulse;  data  for  tests  where  binding  was  suspected  were  dis¬ 
carded. 

3.  Air  drag  forces  are  negligible.  This  was  presumed  to  be  true  because  the  velocities  of 
the  plugs  were  not  high  enough  to  induce  significant  drag  forces  and  because,  for  a 
considerable  portion  of  its  motion,  the  plug  was  surrounded  by  explosion  products 
moving  at  a  higher  velocity. 

The  formulation  represented  by  Equation  (2)  is  identical  with  that  employed  by 

Johnson,  Patterson,  and  Olson/ 


Symbols  are  defined  in  the  List  of  Symbols,  page  37 . 


In  order  to  obtain  sufficient  data  to  map  out  the  variation  of  normally  reflected 
impulse  versus  scaled  distance  it  was  decided  to  perform  tests  at  the  following  values 
of  Z:  • 


m/kg1/3  ft/lb1/3 


0.198 

0.159 

0.119 

0.079 

0.066 

0.060 


0.5 

0.4 

0.3 

0.2 

0.1654* 

0.15 


If  these  were  all  successful  it  is  likely  that  an  attempt  would  have  been  made  to 
measure  I  untU  the  surface  of  the  charge  was  reached;  i.e.,  to  Z  =  0.0525  m/kg1''* 
(0.1323  ft/lb1'3). 

Also,  since  it  was  not  established  that  the  Hopkinson  scaling  law  applies  for  such 
small  Z,  it  was  desired  to  obtain  data  for  more  than  one  charge  weight;  specifically  for 
charges  of  0.227,  0.454,  0.907  kg  (0.5,  1.0,  2.0  lb). 

From  examination  of  standard  deviations  from  earlier  tests  it  was  concluded  that 
three  tests  at  each  parameter-pair  would  yield  results  at  the  95%  confidence  level  while 
the  99%  level  could  be  reached  for  five  replications. 

F.  Testing  Experience 

The  tests  were  conducted  primarily  in  order  of  decreasing  scaled  distance.  Figure  7 
shows  a  representative  photograph  of  the  vertical  scale  with  the  plug  “stopped”  at  two 
positions  by  the  flash  lamps.  Because  the  lamps  are  to  the  right  of  the  camera  the  sha¬ 
dow  of  the  plug  is  to  the  left  of  the  plug.  In  order  to  reduce  the  data  one  must  estimate 
the  vertical  positions  of  the  mass  center;  i.e.,  Xi  and  x2  of  Equation  (2).  This  is  made 
somewhat  more  difficult  by  the  fact  that,  in  spite  of  care  in  centering  the  plug  and  the 
explosive  charge,  most  of  the  plugs  experienced  some  tumbling.  This  tumbling  does  not 
invalidate  the  response  formulation  but  could  increase  the  drag  force.  In  a  few  instances 
the  tumbling  plug  hit  the  scale;  if  this  occurred  above  x2  the  results  of  that  test  were 
discarded. 


This  value  was  chosen  to  correspond  to  the  case  discussed  in  Section  II. 


Figure  7.  Dual  Image  of  Plug  in  Front  of  Scale 

For  the  series  of  tests  with  0.079  kg  (2.0  lbs)  charges  at  Z  =  0.079  m/kg1/3 
(0  *  ft/lb1' )  the  RHA  plugs  were  significantly  deformed  so  that  it  was  necessary  to  use 
th'.  hardened  BEARCAT  plugs.  The  RHA  armor  plate  also  experienced  plastic  flow  and 
'‘i-'.tering,  the  latter  increasing  in  size  and  depth  as  the  tests  continued  at  smaller  Z.  It 
was  found  necessary  to  ream  out  the  cylindrical  hole  in  the  plate  to  its  original  diameter 
before  each  shot.  Another  problem  arose  concerning  the  photographs  of  the  plug 
against  the  scale  background:  fiery  streaks  (burning  explosive  or  metal)  either  blowing 
by  or  following  the  plug  obscured  the  image  of  the  plug  resulting  in  loss  of  data.  Efforts 
to  prevent  this  by  use  of  baffles  were  not  very  successful. 


For  the  smaller  scaled  distances  there  was  concern  that  the  plastic  flow  of  the  plate 
(decrease  of  hole  diameter)  during  the  blast  response  might  result  in  resistance  to  plug 
ejection.  Accordingly,  the  cylindrical  surface  of  the  plugs  was  given  a  thin  coating  of 
manganese  disulfide  powder  in  the  belief  that  its  lubrication  property  would  help  cir¬ 
cumvent  this  possibility. 
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When  the  depth  and  extent  of  the  crater  in  the  armor  plate  became  significant  a 
few  shots  were  fired  in  which  the  crater  was  filled  with  “epoxy  steel.”  This  served  to 
restore  the  original  plate/plug  geometry  at  the  initiation  of  the  shot.  However,  this 
material  was  all  removed  by  the  blast  and  it  was  feared  that  some  of  this  material 
might  intrude  into  the  gap  between  plug  and  plate.  At  the  end  of  this  test  program 
while  attempting  to  obtain  data  for  0.227kg  (0.5  lb)  charges  at  Z  =  0.060  m/kg*' 
(0.15  ft/lb1'*5),  two  shots  were  fired  in  which  hardened  steel  inserts  were  pressed  into 
holes  counterbored  to  a  depth  of  22.2  mm  (0.875  in)in  the  armor  plate.  These  inserts 
were  also  flush  with  the  original  plate  surface  and  the  crater  outside  the  inserts  was 
filled  with  the  epoxy  material.  A  more  positive  retention  method  would  be. required  for 
further  testing  since  the  inserts  were  ejected  during  the  rebound  phase  of  the  plate 
response. 

G.  Impulse  Data 

The  basic  data  obtained  from  the  impulse  plug  tests  are  tabulated  in  y\gpendix  B. 
From  these  data  the  values  of  the  Hopkinson  scaling  parameters  Z  and  I/W"/0  were  cal¬ 
culated  and  are  presented  in  Table  1,  grouped  by  the  nominal  value  of  Z  and,  for  each 
of  these,  by  the  nominal  charge  weight.  Out  of  the  fifty-seven  tests  for  which  complete 
data  were  obtained,  two  were  regarded  with  suspicion  and  their  results  were  rejected  on 
the  basis  of  several  criteria  for  treatment  of  outliers.1  These  tests  are  indicated  by 
asterisks  in  Table  1  and  their  results  were  not  used  in  the  calculations  of  mean  impulses 
or  standard  deviations. 

The  validity  of  the  Hopkinson  scaling  law  was  then  investigated  using  the  statisti¬ 
cal  technique  of  analysis  of  variance.  The  F-statistic  was  employed  at  the  95% 
confidence  level  to  test  the  hypothesis  that  all  the  scaled  impulses  for  each  nominal 
scaled  distance  were  derived  from  the  same  normal  population.  According  to  this  pro¬ 
cedure  one  would  conclude  that  this  hypothesis  is  valid  for 
Z  >  0.159  m/kg1/3  (0.4  ft/lb1/3)  but  should  be  rejected  for  all  smaller  Z  for  which  data 
are  available.  Thus,  it  appears  that  the  Hopkinson  scaling  law  becomes  questionable  as 
the  charge  approaches  the  reflecting  surface.  Such  a  deduction  is  certainly  subject  to 
challenge,  however.  The  charge  weights  have  been  varied  over  only  a  small  range  and 
the  standard  deviations  a  of  the  scaled  impulses  increase  as  Z  decreases.  It  is  unfor¬ 
tunate  that  the  range  of  Z  for  which  a  larger  data  sample  is  required  is  just  the  one 
w'here  the  tests  are  most  difficult  .to  perform.  On  the  other  hand  it  is  noted  that 
for  Z  <  0.159  m/kg1'  (0.4  ft/lb1' the  values  of  I/W1' ^  at  each  Z  exhibit  a  mono¬ 
tone  decrease  as  charge  weight  increases. 

The  values  of  mean  scaled  impulse  (I/W1/3^  are  plotted  against  Z  in  Figure  8,  using 
a  different  symbol  for  each  charge  weight.  Where  they  do  not  conflict  with  other  sym¬ 
bols  “  ±  3  a  bars”  have  been  added  to  data  points  for  the  0.227  kg  (0.5  lb)  charges. 


12 

M.  G.  Natrella,  Engineering  Design  Handbook,  Experimental  Statistics,  Section  4, 
Chapter  17,  AMCP  706-113,  March  1966. 
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V 


.2  (0.5)* 


Test  No. 


0.454  ker  (1  lb)  Pentolite  (approximate! 


.1990  (.5017) 
.1985  (.5004) 
.1980  (.4992) 
.1979  (.4989) 
.1970  (.4982 


Z  =  0.1982  (0.4997) 
I/W1/3  =  7.248  (807.7) 


az  =  0.0006  (0.0014) 
a  =  0.057  (6.4) 


.907kg  (2  lb)  Pentolite  (approximately 


.1982  (.4997) 
.1981  (.4995) 
.1980  (.4992) 
.1989  (.5013) 
.1982  (.4997 


7.17  (799. 
7.07  (788. 
7.29  (813. 
7.37  (821. 
6.94 


Z  =  0.1983  (0.4999)  <rz  —  0.0003  (0.0008) 
1/W1^3  =  7.168  (798.8)  <rx  =  0.171  (19.1) 


TABLE  1.  SCALED  TEST  RESULTS  (CONT’D) 


^nominal  =  0.16  (0.4) 


Test  No. 

Z 

I/W*'0 

0.227  kg  (0.5  lb)  Pentolite  (appr 

oximately) 

98 

.1612  (.4065) 

10.71  (1194.) 

99 

.1588  (.4003) 

10.73  (1196.) 

100 

.1588  (.4002) 

10.98  (1224.) 

101 

.1585  (.3996) 

10.86  (1210.) 

102 

.1585  (.3995) 

10.72  (1194.) 

Z  =  0.1592  (.4012)  <rz  =  0.0012  (0.0030) 

I/W1^3  =  10.800  (1203.6)  =  0.119  (13.3) 

0.454  kg  (1.0  lb)  Pentolite  (approximately) 


38 

.1590  (.4007) 

11.03  (1229.) 

39 

.1588  (.4004) 

10.83(1207.) 

40 

.1588  (.4003) 

10.91  (1216.) 

41 

.1585  (.3995) 

10.90  (1215.) 

42 

.1583  (.3989) 

10.95  (1220.) 

Z  =  0.1587  (0.4000) 

I/W1/3  =  10.925  (1217.5) 


az  —  0.0003  (0.0007) 
(7,  =  0.074  (8.3) 


0.907  kg  (2  lb)  Pentolite  (ap 
.1588  (.4004) 

.1587  (.4001) 

.1586  (.3998) 

.1603  (.4041) 

.1589  (.4005) 


>roximateh 


11.03  (1229.) 
10.71  (1193.) 
10.86  (1211.) 
10.10  (1126.)* 
10.86  (1211.) 


Z  =  0.1588  (0.4002)  <rz  =  0.0001  (0.0003) 
I/W1/3  =  10.865  (1210.8)  or,  =  0.130  (14.5) 


TABLE  1.  SCALED  TEST  RESULTS  (CONT’D) 


^nominal  =  0-12  (0.3) 


Test  No. 

Z 

i/w1/3 

0.227  kg  (0.5  lb)  Pentolite  (approximately) 

103 

.1209  (.3048) 

19.74  (2200.) 

104 

.1192  (.3006) 

19.97  (2226.) 

107 

.1183  (.2983) 

19.82  (2208.) 

109 

.1187  (.2992) 

19.54  (2178.) 

110 

.1186  (.2989) 

19.12  (2131.) 

Z  =  0.1191  (0.3003)  <rz  =  0.0010  (0.0026) 

I/W1/3  _  10  640  (2188.6)  ax  —  0.329  (36.6) 


0.454  kg  (1.0  lb)  Pentolite  (approximately) 


43 

.1196  (.3015) 

18.86  (2101.) 

44 

.1196  (.3014) 

18.40  (2050.) 

45 

.1192  (.3005) 

18.51  (2063.) 

52 

.1190  (.3000) 

17.35  (1934.) 

54 

.1195  (.3012) 

18.92  (2109.) 

56 

•1192  (.3005) 

19.04  (2122.) 

Z  =  0.1194  (0.3010)  <rz  =  0.0092  (0.0005) 

I/W1/3  =  18.745  (2089.0)  <rr  =  0.276  (30.8) 


0.907  kg  (2.0  lb)  Pentolite  (approximately) 


76 

.1200  (.3026) 

18.04  (2011.) 

77 

.1194  (.3009) 

18.94  (2110.) 

78 

.1190  (.3001) 

18.78  (2092.) 

79 

.1190  (.2999) 

18.78  (2093.) 

80 

.1189  (.2997) 

18.83  (2098.) 

Z  =  0.1193  (0.3006)  trz  =  0.0005  (0.0012) 

I/W1/3  =  18.673  (2080.9)  <r,  =  0.358  (39.8) 


TABLE  1.  SCALED  TEST  RESULTS  (CONT’D) 


^nominal —  (®‘2) 


- 

Z 

i/w1/3 

0.227  kg  (0.5  lb)  Pentolite  (approximately) 

.0790  (.1992) 

419  (4670.) 

.0790  (.1991) 

41.8  (4660.) 

.0794  (.2003) 

45;2  (5040.) 

.0794  (.2001) 

43.6(4860.) 

.0794  (.2001) 

45.7  (5090.) 

Z  =  0.0792  (0.1997)  <rz  =  6.0002  (0.0006) 
I/W1/3  =  43.65  (4865.)  <rx  =  1.81  (202.) 


0.454  kg  (1.0  lb)  Pentolite  (approximately) 

.0796  (.2006) 

41.0  (4570.) 

.0795  (.2003) 

42.3(4710.) 

.0793  (.1999)- 

42.6  (4740.) 

.0799  (.2014) 

42.2  (4710.) 

.0796  (.2008) 

43:2(4810.) 

Z  =  0.0796  (0.2006)  <rz  =  0.0002  (0.0006) 

I/W1/3  =  42.26  (4709.)  Oj  =  0.78  (86.5) 


0.907  kg  (2.0  lb)  Pentolite  (approximately) 

.0794  (.2001) 

40.4  (4500.) 

.0793  (.1999) 

40,3  (4500.) 

.0793  (.1998) 

39.0  (4340.) 

Z  =  0.0793  (0.1999)  <rz  =  0.0006  (0.0002) 

I/W1''3  =  39.90  (4446.)  ffj  =  0.80  (89.0) 


TABLE  1.  SCALED  TEST  RESULTS  (CONT’D) 


Curves  have  been  drawn  through  the  data  points  for  the  smallest  and  largest  charges 
tested.  These  curves  are  seen  to  coalesce  for  the  Z  range  where  the  Hopkinson  law  was 
validated  by  the  statistical  analysis  and  to  fair  in  nicely  with  the  data  of  Johnson, 
Patterson,  and  Olson  ior  larger  Z. 


i 

Z  (ft/lb3) 


Figure  8.  Scaled  Impulse  Versus  Scaled  Distance 


One  could  certainly  draw  a  single  “best-fit”  curve  for  all  the  data  of  this  report; 
however,  this  could  be  misleading  if  significantly  different  charge  weights  were  to  be 
used.  While  the  cause  of  this  apparent  dependence  of  scaled  impulse  on  charge  weight 
for  small  Z  has  not  been  identified  it  should  be  noted  that  geometric  scaling  was  not 
preserved  in  these  experiments  in  that  no  attempt  was  made  to  vary  the  diameter  of  the 
plug  as  Z  was  varied.  Figure  9  compares  the  geometries  of  the  charge/plug 
configurations  at  the  extremes  of  charge,  weights  and  smallest  Z  tested. 


(2  lb)  CHARGE 


(1/2  lb)  CHARGE 


WMv. 

it 


Figure  9.  Comparison  of  Geometries  at  Z  *=  0.060  m/kg  '  (0.15  ft/lb  '  ) 


H.  Pressure  Measurements 

During  this  investigation  three  exploratory  tests  were  performed  in  which  reflected 
overpressures  were  recorded.  Two  PCB  Model  No.  109A  piezoelectric  transducers  were 
mounted  flush  with  the  upper  surface  of  the  armor  plate  and  152.4  mm  (6  in)  apart. 
The  spherical  charge  was  positioned  directly  above  one  of  these  gages.  Figure  10  and  11 
show  the  recorded  pressures  from  these  gages  for  two  firing  of  0.454  kg  (1  lb)  charges  at 
a  height-of-burst  corresponding  to  Z  =  0.2  m/kg1' **  (0.5  ft/lb1'  ).  While  some  evi¬ 
dence  of  ringing  of  the  armor  plate  may  be  present  in  these  records,  they  serve  to  pro¬ 
vide  insight  into  the  complexities  of  blast  wave  reflections  at  small  Z.  Figure  10  ivus- 
trates  the  variability  of  wave  forms  for  normal  reflection  which  nmy  occur  for  nominally 
identical  firings.  For  this  case  the  tables  of  Kingerv  and  Pannilr  gave  a  peak  reflected 
overpressure  of  168.68  MPa  (24465  psi)  while  Jack1,1  has  reported  the  value  186.4  MPA 


5W.  H.  Jack,  Jr.,  “Measurements  of  Normally  Reflected  Shock  Waves  from  Explosive 
Charges,”  Ballistic  Research  Laboratories  Memorandum  Report  No.  1499,  July  1963. 


(27040  psi)  as  the  average  of  three  firings  using  0.057  kg  (1/8  lb)  charges.  In  Figure  11 
the  reflected  overpressure  records  for  the  same  two  firings  are  seen  to  be  in  much  closer 
agreement.  The  location  of  this  gage  with  respect  to  “ground  zero”  is  just  outside  of 
the  regular  reflection  zone. 

For  the  same  charge  weight  a  single  firing  was  made  at  Z  =  0.12  m/kg1/3 
(0.3  ft/lb1'3);  the  blast  overpressure  recorded  at  the. same  two  gage  locations  are  shown 
in  Figures  12  and  13,  respectively.  The  very  high  pressures  shown  in  Figure  12  resulted 
in  damage  to  this  gage;  it  is  not  clear  whether  the  true  peak  reflected  pressure  was 
reached  before  this  occurred.  According  to  Kingery  and  Pannill'  one  would  eraect  a 
peak  reflected  pressure  of  only  348.0  MPa  (50,476  psi).  However,  Kulesz  et  al  have 
reported  an  experimental  value  of  830  MPa  (120,000  psi)  for  this  case.  Interpolating 
data  from  this  source  one  would  estimate  a  peak  pressure  of  210  MPa  (30,000  psi)  for 
the  location  of  the  gage  output  shown  in  Figure  13,  which  is  clearly  not  the  result 
obtained.  Owing  to  the  decreased  height-of-burst  this  location  is  a  bit  further  into  the 
Mach  reflection  region  than  was  the  case  for  Figure  11. 

The  blast  pressure  measurements  reported  in  Figures  10-13  are  too  few  in  number 
to  have  any  statistical  significance  but  do  indicate  that  uncertainties  associated  with 
blast  load  predictions  increase  rapidly  as  Z  is  reduced.  More  importantly,  it  should  be 
noted  that  none  of  these  blast  pulses  have  the  sharp  rise  followed  by  exponential  decay 
that  is  typical  at  larger  scaled  distances.  These  anomalous  wave  forms  pose  serious 
problems  in  the  evaluation  of  such  blast  parameters  as  positive  phase  duration  and  posi¬ 
tive  impulse. 


The  determination  of  the  reflected  positive  duration  was  a  prime  objective  of  these 
pressure  measurements  since  this  is  required  for  estimation  of  plug  motion  during  the 
blast  pulse.  The  displacement  at  the  end  of  the  blast  pulse  is  given  by  Equation  (A-6) 
where,  for  short  pulses,  the  first  term  is  negligible  in  comparison  to  the  second.  C  is  a 
non-dimensional  pulse  shape  factor  which,  for  the  shapes  considered  in  Appendix  A,  is  in 

the  range  —  <  C  <  1.  Thus,  for  a  specified  velocity  x0,  the  displacement  x(T)  is 
2 

directly  proportional  to  T.  Although  data  on  positive  duration  for  small  Z  are  rather 
sparse  it  appears14  that  this  quantity  is  nearly  constant  over  the  range  0.08  <  Z  .<  0.20 
m/kg y  3  (0.15  <  Z  <  0.50  ft/lb1/3)  and  is  of  the  order  of  0.13  ms/kg1'  (0.1 
ms/lb1'  ),  consistent  with  the  positive  durations  seen  in  Figures  10-13.  On  this  basis  it 
is  felt  that  the  plug  motion  during  the  blast  pulse  may  have  been  somewhat  greater 
than  anticipated  for  tests  at  the  smallest  scaled  distances. 


14C.  N.  Kingery  and  G.  Bulmash,  “Airblast  Parameters  from  TNT  Spherical  Air  Burst 
and  Hemispherical  Surface  Burst,”  Ballistic  Research  Laboratory  Technical  Report 
ARBRL-TR-02555,  April  1984. 


Z  =  0.12  m/kg3  (0.3  ff/lb3) 
R  =152.4  mm  (6 in) 
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IV.  CONCLUDING  REMARKS 


Data  on  normally  reflected  impulse  from  spherical  Pentolite  charges  have  been 
presented  for  the  range  of  small  scaled,  distances  not  adequately  treated  in  the  literature. 
Although  the  experimental  measurements  are  believed  to  be  normally  distributed  the 
effects  omitted  in  the  interpretive  analysis:  mechanical  friction,  air  drag,  and  pressure 
relieving  motion  of  the  plug  during  the  positive  phase  of  the  blast  pulse  all  would  tend 
to  reduce  the  plug  velocity.  Consequently,  the  reported  values  should  represent  a  (hope¬ 
fully  close)  lower  bound  to  the  normally  reflected  impulse  which  an  infinite  rigid  plane 
would  experience. 

There  is  tentative  evidence  that  the  Hopkinson  scaling  law  (and  consequentially, 
Sachs’  scaling)  may  become  inapplicable  for  Z  <  0.16  m/kg1'  (0.4  ft/lb1'  ).  Further 
evaluation  of  this  matter  using  a  greater  variation  of  charge  weights  is  needed. 

The  experimental  difficulties  encountered  in  this  program  were  described,  as  were 
the  corrective  measures  which  were  undertaken.  For  future  testing  a  reliable  method  for 
measuring  the  transient  velocity  of  the  plug  should  be  established  as  well  as  a  means  of 
determining  the  pulse  duration.  Provision  for  simultaneous  recording  of  surface  pres¬ 
sures  as  a  function  of  distance  R  should  be  made;  these  are  needed  for  structural 
response  calculations. 

With  respect  to  the  comparison  of  impulses  shown  in  Figure  2,  the  experimental 
value  from  Test  No.  89  may  be  included  at  R/H  —  0.  From  Table  R-2  one  finds 
I  =  51.3  kpa-s  (7440  psi-ms),  which  is  somewhat  less  than  the  empirical  values  shown 
in  Figure  2.  Although  there  is  only  a  single  data  point  for  this  case  it  may  be  seen  from 
Figure  8  that  it  is  consistent  with  those  for  neighboring  values  of  Z.  Therefore  it  may 
be  concluded  that  the  normally  reflected  impulse  derived  from  the  HULL  hydrocode  is 
not  more  than  60%  of  the  best  available  experimental  value.  It  is  hoped  that 
refinements  to  the  modeling  of  the  detonation  process  and  of  the  treatment  of  explosive 
products  can  be  devised  which  will  bring  the  hydrocode  predictions  into  line  with  the 
experimental  evidence.  This  is  sorely  needed  since  computational  hydrodynamics  can 
provide  detail  regarding  spatial  distribution  of  pressure  and  impulse  which  is  infeasible 
to  obtain  experimentally. 

ACKNOWLEDGMENTS 

Numerous  individuals  provided  assistance  during  the  conduct  of  this  investigation' 
and  the  authors  wish  to  express  their  thanks  to  all  who  participated.  R.  Lottero  pro¬ 
vided  the  HULL  code  solution  referenced  herein.  M.  Saccucci  helped  with  the  statistical 
design  of  the  experiments.  W.  T.  Robinson  assisted  with  performance  of  the  experi¬ 
ments  and  with  collection  of  information  to  be  included  in  the  report.  R.  Schumacher 
provided  support  for  data  reduction  and  plotting.  C.  N.  Kingery  provided  extrapolated 
data  from  experiments  performed  by  Southwest  Research  Institute  and,  during  the 


33 


15 

review  of  this  report,  called  attention  to  a  relevant  paper  by  Baker  which  is  taken 
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with  the  principal  author. 
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LIST  OF  SYMBOLS 


A  area  of  (circular)  cross  section  of  plug 
C  pulse  shape  factor 

H  height-of-burst  (to  center  of  sphere) 

I  impulse  per  unit  area  , - 

K  coefficient  in  Equation  (C-2),  =  ^[2Y/Mg/47Tg 
Mg  mass  of  explosive  charge 

Ma  mass  of  air  engulfed  by  shock  wave 

R  horizontal  radial  distance  to  “ground  zero” 

T  duration  of  positive  phase  of  blast  pulse 
W  weight  of  spherical  explosive  charge 

Y  energy  released  by  Mg 

Z  Hopkinson  scaled  distance  ( =  H/W1/3) 
g  acceleration  of  gravity 

m  mass  of  plug 

p  transient  reflected  pressure 

pQ  peak  reflected  pressure 

t  time 

x  vertical  position  of  mass  center  of  plug  relative 

to  its  initial  position 

XQ  velocity  imparted  to  plug  by  blast  alone  (=  Al/m) 
1 3  decay  coefficient  of  the  Friedlander  pulse 

7  ratio  of  specific  heats  (=  cp  /  cv) 
a  standard  deviation 

T  dummy  time  variable 


APPENDIX  A 

ANALYSIS  OF  IMPULSE  PLUG  RESPONSE 


If  the  effect  of  air  resistance  is  neglected,  the  motion  of  the  plug  in  the  presence  of 
the  gravitational  field  is  described  by  the  differential  equation 

mx  =  A  p(t)  +  mg  ,  (A- 1) 

subject  to  the  initial  conditions 

x(0)  =  0  ,  X  (0)  =  0  (A-2) 

Equation  (A-l)  may  be  successively  integrated  to  obtain  expressions  for  the  velocity  and 
the  displacement: 

x(t)  =  gt  -f  — / 1  p(r)dr  (A-3) 

m  o 

x(t)  =  ygt2  +  -^-jf  *  p(r)  (t  -  r)  dr  (A-4) 

In  order  to  proceed  further  one  must  be  more  specific  regarding  the  form  of  p(t).  Three 
possible  pulse  shapes  were  considered  in  this  investigation: 

(1)  a  linearly  decaying  (triangular)  pulse,  which  had  been  employed 
by  Johnson,  Patterson,  and  Olson; 

(2)  the  Friedlander  function,  which  is  often  used  to  represent  the 
classical  decay  of  a  blast  wave;  and 

(3)  a  half-sine  pulse,  which  was  included  iii  view  of  the  non-elassieal 
pressure  records  observed  experimentally  for  small  Z. 


Figure  A-I  presents  sketches  of  the  pulses  and  the  associated  analytical  formulations  for 
the  transient  pressures. 

For  any  blast  pulse  where  the  response  of  the^plug  is  described  by  Equation  (A-l) 
the  velocity  at  the  end  of  the  pulse  is: 


x(T)  =  gT  + 


m 


(A- 5) 


AI 


Let  XQ  =  - be  the  velocity  imparted  to  the  plug  by  the  blast  pulse  in  the  absence  of 

m 

a  gravitational  field.  Then,  the  displacement  of  the  plug  at  the  end  of  the  pulse  is  given 
by: 


x  (T)  =  {gT2  +  CT  i0 


(A- 6) 


where  the  dependence  on  pulse  shape  is  incorporated  in  the  factor  C  which  has  the 

2  2-2/5  +  f-  2e^  1 

values  —  ,  — S' In - ~ —  >  7T  >  respectively,  for  the  three  pulses  presented 

O 

in  Figure  A- 
are  given  by: 


i.  it 


-  l  +  /?)  ’  2 


e  velocities  and  displacements  of  the  plug  after  the  end  of  the  pulse 


\2r-jW 


A.*  »1»V 


cl 


X  =  gt  +  xc 


x  =  -|gt2  +  x0t  +  (C-l)x„T 


t  >  T 


(A-7) 


(A-8) 


Equation  (A-8)  provides  expressions  for  plug  locations  Xj  ,  x2  at  times  tj  ,  t2, 
where  t2  >  tj  >  T.  Taking  the  difference  of  these  expressions  one  obtains: 


X0  -  X 


to  -  t 


- *7  (fc2  +  tl) 


(A-9) 


This  relation,  while  valid,  has  the  disadvantage  that  the  times  fcj  and  t2  from  the 

arrival  of  the  blast  pulse  at  the  plug  must  be  known.  Since  it  is  easier  experimentally  to 
measure  t2  -  fcj,  write: 


*1  =  gtj  +  x0 

and  substitute  for  xQ  from  (A-9)  to  obtain: 


(A- 10) 


X,  = 


X2-  x 


to  —  t 


"T  (^2  -  t.l) 


(A-ll) 


Time  can  be  eliminated  between  Equations  (A-7)  and  (A-8)  (written  for  t 
obtain  a  quadratic  expression  in  XQ  which  yields: 

*o  =  V^i2  "  2gxj  +  (1  -  C)2  g2  T2  -  (1  -  C)  gT 


=  tj)  to 


(A-12) 


In  order  to  employ  Equation  (A-12).  one  needs  the  value  of  the  pulse  duration  T  (which 
was  not  measured  in  th|^lug  experiments).  Data  on  normally  reflected  positive  dura¬ 
tion  published  by  Jack  indicate  a  monotonic  decrease  of  duration  with  Z  but  no 
values  were  eiven  for  Z  fl.4  m/kcr*/3  /%  _ „i _ j  . 


tion  published  by  Jack  indicate  a  monotonic  decrease  of  duration  with  Z  but  no 
values  were  given  for  Z  <  0.4  m/kg1/3  (1  ft/lb1/3).  The  durations,  observed  in  the 
three  pressure  measurement  tests  reported  herein  were  all  a  small  fraction  of  a  mil¬ 
lisecond.  A  study  was  conducted  using  representative  data  from  this  report  in  Equation 
(A-12),  taking  T  to  be  both  0  and  1  ms.  The  difference  in  the  velocity  was  found  to  be 
of  the  order  of  0.01%  or  less  so  it  was  concluded  that  the  terms  involving  T  could  be 
omitted  from  Equation  (A-12)  for  the  present  application.  On  this  basis  the  value  of 
impulsive  velocity  is  given  by: 

X6  =  V*!2  -  2gXj  (A- 13) 


with  the  value  of  Xj  obtained  from  Equation  (A-ll).  These  results  provide  the  basis  for 
Equation  (2)  of  this  report. 
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APPENDIX  B 


EXPERIMENTAL  MEASUREMENTS  AND  UNSCALED  RESULTS 


APPENDIX  B 

EXPERIMENTAL  MEASUREMENTS  AND  UNSCALED  RESULTS 

Table  B-l  is  a  listing  of  the  basic  data  obtained  in  each  test,  given  in  the  units 
actually  employed  in  the  measurements.  Results  of  some  intermediate  calculations  are 
listed  in  Table  B-2.  Tabulated  in  both  metric  S.I.  and  English  units  are  the  plug  velo¬ 
city  at  the  end  of  the  pulse  and  the  corresponding  impulse  (per  unit  area). 


TABLE  B-l.  SUMMARY  OF  TEST  DATA 


Test  No. 

Charge 

Weight 

W 

grams 

H 

in. 

Plug 

Weight 

grains 

Plug 

Diameter 

in. 

Xl 

in. 

x2 

in. 

ms. 

h 

ms 

33 

479.7 

6.134 

384.5 

.4985 

13.8 

31.3 

12.58 

28.04 

34 

483.5 

6.134 

384:5 

.4985 

13.9 

31.6 

12.52 

27.90 

35 

487.0 

6.134 

384.5 

.4985 

13.6 

31.3 

12.32 

27.72 

36 

488.0 

6.134 

384.5 

.4985 

13.7 

31.2 

12.40 

27.78 

37 

490.0 

6.134 

384.5 

.4985 

13.7 

31.2 

12.60 

28.02 

38 

484.4 

4.915 

384.5 

.4985 

14.2 

32.3 

8.66 

19.13 

39 

485.5 

4.915 

384.5 

.4985 

143 

32.2 

8.60 

19.14 

40 

485.9 

4.915 

384.5 

.4985 

14.2 

’  32.3  • 

8.69 

19.26 

41 

488.9 

4.915 

384.5 

.4985 

14.3 

-  32.4 

8.65 

19.21 

42 

490.9 

4.915 

384.5 

.4985 

14.2 

32.4 

8.67 

19.23 

43 

480.0 

3.687 

384.5 

.4985 

13.6 

33.8 

4.86 

11.73 

44 

480.4 

3.687 

384.5 

.4985 

13.9 

33.7 

•4.98 

11.88 

45 

484;8 

3.687 

384.5 

.4985 

14.1 

34.0  . 

4.97 

1184 

52 

487.3 

3.687 

388. 

.4995 

148 

30.9 

5.64 

11.59 

.  54 

481.7 

3.687 

388. 

.4995 

1519 

33.3 

5.7Q 

11.62 

56 

484.9 

3.687 

388. 

.4995 

15.8 

34.0 

5.69 

11.83 

65 

896.5 

7.5252 

389. 

.4995 

16.3 

35.5 

12.38 

26.38 

66 

897.7 

7.5252 

389. 

.4995 

13.5 

30.8 

10.48 

23.26 

67 

899.1 

7.5252 

389. 

.4995 

13.8 

31.5 

10.38 

23.06 

68 

887.7 

7.5252 

389. 

.4995 

14.1 

31.8 

10.62 

23.22 

69 

896.6 

7.5252 

389. 

.4995 

13.3 

30.0 

10.64 

23.22 

70 

893.8 

6.024 

389. 

.4995 

13.1 

31.8 

6.64 

15.54 

71 

896.2 

6.024 

389. 

.4995 

13.0 

31.3 

6.68 

15.64 

72 

898.1 

6.024 

389. 

.4995 

13.2 

31.8 

6.70 

15.67 

74 

869.9 

6.024 

389. 

.4995 

12.4 

29.4 

6.66 

15.57 

75 

893.3 

6.024 

389. 

.4995 

13.0 

31.4 

6.66 

15.55 

47 
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TABLE  B-l.  SUMMARY  OF  TEST  DATA  (CONT’D) 


Test  No. 

Charge 

Weight 

W 

crams 

H 

in. 

Plug 

Weight 

strains 

Plug 

Diameter 

Xl 

in. 

*2 

in. 

76 

872.2 

4.515 

389. 

.4995 

15.1 

29.1 

77 

886.4 

4.515 

389. 

.4995 

16.0 

30.7 

78 

894.2 

4.515 

389. 

.4995 

16.3 

30.9 

79 

895.6 

4.515 

389. 

.4995 

16.1 

30.8 

80 

897.4 

4.515 

389. 

.4995 

16.0 

30.8 

83 

893.0 

3.010 

380. 

.4968 

35.3 

52.25 

84 

8P6.1 

3.010 

380. 

.4968 

36.0 

52.8  : 

87 

897.5 

3.010 

386. 

.4997 

36.5 

52.3 

89 

895.6 

2.489 

387. 

.4985 

33.2 

47.4 

92 

485.0 

2.461 

384. 

.4978 

36.0 

51.7 

93 

487.0 

2.461 

384. 

.4978 

34.1 

53.5 

95 

489.9 

2.461 

384. 

.4978 

30.5 

50.2 

96 

479.1 

2.461 

384. 

.4978 

30.0 

49.6 

97 

483.5 

2.461 

384. 

.4978 

30.8 

50.9 

98 

227.8 

3.877 

384. 

.498 

28.7 

43.7 

99 

238.5 

3.877 

384. 

.498 

29.1 

44.6 

100 

238.6 

3.877 

384. 

.498 

30.4 

46.1 

101 

239.8 

3.877 

384. 

.498 

29.5 

45.1 

102 

240.0 

3.877 

384. 

.498 

28.7 

44.1 

103 

228.0 

2.908 

384. 

.498 

31.0 

46.7 

104 

237.7 

2.908 

384. 

.498 

31.3 

47.6 

107 

243.3 

2.908 

384. 

.498 

31.5 

47.7 

109 

241.1 

2.908 

384. 

.498 

31.4 

17.4 

110 

241.7 

2.908 

384. 

.498 

31.5 

47.1 

111 

241.9 

1.938 

384. 

.498 

34.5 

51.8 

112 

242.2 

1.938 

384. 

.498 

29.8 

47.0 

113 

227.4 

1.909 

384. 

.498 

30.5 

48.6 

114 

227.9 

1.909 

384. 

.498 

29.9 

47.4 

115 

228.0 

1.909 

384. 

.498 

31.8 

50.2 

116 

238.7 

1.454 

383. 

.497 

32.8 

48.1 

125 

242.3 

1.454 

385. 

.498 

30.2 

44.9 

4.70 

4.74 

4.765 

4.70 

4.72 

4.63 

4.64 
4.735 
3.05 
5.715 
5.24 
4.705 
4.74 

4.73 
23.02 
23.12 
23.04 
23.16 
23.16 

13.47 
13.46 

13.48 
13.65 
13.64 

6.85 

5.865 

5.80 

5.82 

5.84 

3.955 

3.95 


8.81 

8.83 

8.85 

8.81 

8.845 
6.81 
6.80 

6.845 
4.45 
8.165 
8.175 
7.66 
7.725 
7.715 

34.46 

34.74 

34.54 

34.7 

34.7 

20.0 

20.07 

20.05 

20.25 

20.21 

10.18 

9.185 

9.095 

9.12 

9.155 

5.71 

5.70 
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TABLE  B-2.  PLUG  VELOCITIES  AND  UNSCALED  IMPULSES 


Test  No. 


33 

34 

35 

36 

37 

38 

39 

40 

41 

42 

43 

44 

45 
52 
54 
56 

65 

66 

67 

68 

69 

70 

71 

72 

74 

75 

76 

77 

;s 

79 

80 

83 

84 
87 
89 


5.65 
5.75 
5.74 

5.68 
5.67 

8.66 

8.51 

8.58 

8.59 

8.64 
14.76 
14.41 
14.54 
13.65 
14.93 
14.96 

6.91 

6.82 

7.04 

7.08 

6.69 
10.62 
10.32 
10.48 

9.64 
10.46 
17.24 
18.19 
18.09 
18.10 
18.16 
38.9 
38.9 
37.6 
51.3 


820. 

833. 

832. 

824. 

822. 

1257. 

1234. 

1245. 

1246. 
1253. 
2140. 
2090. 
2110. 
1980. 
2150. 
2170. 
1002. 

990. 

1021. 

1027. 

971. 

1540. 

1497. 

1520. 

1399. 

1517. 

2500. 

2640. 

2620. 

2620. 

2630. 

5640. 

5640. 

5450. 

7440. 


TABLE  B-2.  PLUG  VELOCITIES  AND  UNSCALED  IMPULSES  (CONT’D) 


Test  No. 

y. 

m/s 

0 

in/ms 

kpa-s 

I 

psims 

92 

162.7 

6.41 

32.2 

4680. 

93 

167.8 

6.61 

33.3 

4820. 

95 

169.3 

6.66 

33.5 

4870. 

96 

166.7 

6.56 

33.0 

4790. 

97 

171.0 

6.73 

33.9 

4910. 

98 

33.0 

1.300 

6.54 

949. 

99 

33.6 

1.323 

6.66 

965. 

100 

34.4 

1.354 

6.81 

988. 

101 

34.1 

1.341 

6.75 

978. 

102 

33.6 

1.324 

6.66 

966. 

103 

60.9 

2.40 

12.06 

1749. 

104 

62.5 

2.46 

12.37 

1795. 

107 

62.5 

2.46 

12.37 

1794. 

109 

61.4 

2.42 

12.16 

1764. 

110 

60.1 

2.37 

11.91 

1728. 

111 

131.9 

5.19 

26.1 

3790. 

112 

131.5 

5.18 

26.0 

3780. 

113 

139.5 

5.49 

27.6 

4010. 

114 

134.6 

5.30 

26.7 

3870. 

115 

140.9 

5.55 

27.9 

4050. 

116 

221. 

8.72 

43.9 

6370. 

125 

213. 

8.40 

42.4 

6140. 
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.  APPENDIX  C 

CORRELATIONS  WITH  PREDICTIONS  EMPLOYING  BAKER’S  MODEL 


W.  E.  Baker13  recommends  the  semi-empirical  formula 

I  =  ^  2(Me  +  MA)  Y  /  4  7T  H2  (C-l) 

for  prediction  of  the  normally  reflected  impulse  imparted  to  a  rigid  wall.  He  also  shows 
that  for  Ma  <  <  Mg.  Equation  (C-l)  is  consistent  with  Hopkinson’s  scaling  law  and 
that  a  reasonable  upper  limit  on  Z  for  neglect  of  MA  is  0.269  m/kg1' 3  (0.679  ft/lb1'3), 
if  ma  can  be  neglected  Equation  (C-l)  can  be  readily  manipulated  into  the  form 

I  /  W1/3  =  K  /  Z2  (C-2) 

where  K  is  a  dimensional  coefficient  which  can  be  adjusted  to  match  data  for  impulse 
delivered  by  specific  explosives.  Using  the  data  of  this  report  for  50/50  Pentolite  at  Z 
=  0.2  m/kg1' °  (0.5  ft/lb  '  )  one  obtains: 

0.2833  /  Z2  in  S.  I.  units 

200.65  /  Z2  in  English  units 

This  simple  inverse  square  dependence  on  Z  has  been  employed  to  predict  values  of 
scaled  impulse  coi  responding  to  the  mean  experimental  values  shown  in  Table  1.  The 
values  of  I/W1'"  obtained  from  Equation  (C-3)  and  the  percent  differences  from  the 
corresponding  experimental  data  are  listed  in  Table  C-l.  It  may  be  seen  that  Baker’s 
procedure  provides  values  of  reflected  impulse  which  are  quite  satisfactory  except  for  the 
very  small  Z  range  where  it  has  already  been  concluded  that  Hopkinson  scaling  may 
become  inapplicable. 


TABLE  C-l.  SCALED  IMPULSES  DERIVED  FROM  EQUATION  (C-3) 


z 

ft/lb1/3 

Nominal 
Charge  Wt. 

lb 

i/w1/3 

lcDa*s/ke^3  psi°ms/lb*/3 

%Difference 
from  Exp.  Value 

0.1982 

EH 

1.0 

7.212 

803.6 

-0.5 

0.1983 

mm 

2.0 

7.204 

802.9 

0.5 

0.1592 

0.4012 

0.5 

11.178 

1246.6 

3.5 

0.1587 

0.4000 

1.0 

11.248 

1254.1 

3.0 

0.1588 

0.4002 

2.0 

11.234 

1252.8 

3.4 

0.1191 

0.3003 

0.5 

19.972 

2225.0 

1.7 

0.1194 

0.3010 

1.0 

19.872 

2214.7 

6.0 

0.1193 

0.3006 

2.0 

19.905 

2220.6 

6.7 

0.0792 

0.1997 

0.5 

45.16 

5031. 

3.4 

0.0796 

0.2006 

1.0 

44.71 

4986. 

5.8 

0.0793 

0.1999 

2.0 

45.05 

5021, 

12.9 

0.0656 

0.1653 

2.0 

65.83 

7343. 

23.8 

0.0594 

0.1497 

0.5 

80.29 

8954. 
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ADTC  (Tech  Lib) 

Eglin  AFB,  FL  32542-5030 

AFATL 

(DRLRV,  J.  R.  Rutland) 
(DLYV) 

Eglin  AFB,  FL  32542-5300 

USAFTAWC  (3A) 

Eglin  AFB,  FL  32542-5030 

TAC 

Langley  AFB,  VA  23365 

AFAPL  (3FH) 

Wright-Patterson  AFB, 

OH  45433 

AFFDL  (FES) 

Wright-Patterson  AF3, 

OH  45433 

AFFDL  (FDTR,  Dr.  JaniiO 
Wright-Patterson  AFB, 
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1  AFFDL  (F3E,  Dr.  Bader) 
Wright-Patterson  AFB, 

OU  45433 

1  AFLC  (MCNEA) 

Wright-Pattersoh  AFB, 

ON  45433 

1  AFLC  (PPSC) 

Wright-Patterson  AFB, 

OH  45433 

1  AFLC  (MMhiO 

Wright-P.atterson  AFB, 

OH  45433 

1  AFLC  (HEA) 

Nright-Patterson  AFB, 

OH  45433 

1  AFML  (LLN,  Dr.  Nicholas) 

Wright-Patterson  AFB,  OH  45433 


2  Director 

SanJia  National  Laboratories 
ATTN:  Doc  Cont  for  3141 

SanJia  Report  Collection 
L.  J.  Vortman 
Albuquerque,  NM  37115 

1  Director 

NASA  Scientific  and  Tech  Info 
Facility 
P.O.  Box  3757 
BaltWash  international 
Airport,  HD  21240 

1  Wilfred  Baker  Engineering 
213  £.  Edgewood  Place 
San  Antonio,  TX  73209 

1  Aerospace  Corporation 

ATTN:  Tech  Info  Services 

P.O.  Box  92957 

Los  Angeles,  CA  90009 


1  AFIT  (Lib  Bldg.  540,  Area  3)  2  Battelle  Colunbus  Lab 

Wright-Patterson  AFB  OH  45433  ATTN:  Dr.  L.  E.  Hulbert 

Tech.  Library 

1  FTD  (TD/BTA/Lio)  505  Xing  Ave 

Wright-Patterson  AFB,  OH  45433  Colunbus,  OH  43201 
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Strategic  Air  Command 
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Lawerence  Livermore  Lao 
Tech  Info  Div 
P.O.  Box  303 
Livermore,  CA  94550 
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Los  Alamos  Scientific  Lab 
ATTN:  Dae  Control  for  Reports 
Library 

C.  L.  Mader 
P.O.  Box  1653 
Los  Alanos,  H.1  37544 


1  Director 

Lewis  Directorate  4  Development 
Lab 

Lewis  Research  Center 
(Hail  Stop  77-5) 

Cleveland,  OH  44135 

1  Bell  Helicopter  Textron 
ATTN:  Hr.  Nile  Fischer 
P.O.  Box  432 
Fort  Worth,  TX  75101 

1  The  3DM  Corporation 
ATTN:  Richard  Hensley 
P.O.  Box  9274 
Albuquerque  International 
Albuquerque,  NM  37119 
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1  Ine  Boeing  Company 
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P.0.  3ox  3999 
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1  Tne  Boeing  Company 
Vertol  Division 
ATTN:  Dave  Harding 
P.0.  Box  16353 
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1  Goodyear  Aerospace  -Corp 
ATTN:  R.  M.  3rown,  Bldg  1 

Shelter  Engineering 
Litchfield  Park,  AZ  35340 

2  Denver  Research  Institute 
University  of  Denver 
ATTN:  Mr.  J.  Vfisotski 

Technical  Library 
2390  South  University  31vd 
Denver,  CO  33210 

1  Fairchild  Industries 
Fairchild  Republic  Division 
ATTN:  D.  Nat son 
Farmingdale  LI,  NY  11735 

2  Falcon  Research  and  Development 

Company 

ATTN:  L.  Mahood 

A.  Stein 
109  Inverness  Drive 
Englewood,  CO  30112 

2  Xaman  AviDyne 

ATTN:  Dr.  N.  P.  Hobbs 

Mr.  3.  Criscione 
33  Second  Avenue 
Northwest  Industrial  park 
Burlington,  MA  01303 
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1  Kaman-TEMPO 
ATTN:  DASIAC 

P.Q.  Drawer  QQ,  31 6 State  St 
Santa  Barbara,  CA  93102 

1  Kaman  Tempo 

ATTN:  E.  Bryant 
H  5.  Shamrock  Road 
Suite  UL-1 
Bel  Air  ,  MD  21 01 4  - 

1  Hughes.  Helicopters 
31dg  5  M/S  C303 

ATTN:  Security  Officer 
Centinela  and  Teale  Sts 
Culver  City,  CA  90230 

2  University  of  Delaware 
Dept  of  Mech  i  Aerospace 

Engineering 

ATTN:  Prof.  J.  R.  Vinson 
Prof,  M.,  Taya 
Newark,  D£  19711 

1  Florida  Atlantic  U’niv 
Dept  of  Ocean  Eng 
ATTN:  Prof.  K.  Stevens 
Boca  Raton,  FL  33432 

1  J.G.  Eng  Research  Associates 
3331  Menlo  Drive 
Baltimore,  MD  21215 

1  Mass  Inst  Of  Tech 

Aeroelastic  1  Structures 
Research  Laboratory 
ATTN:  Dr.  E.  A.  Ni  truer 
Cambridge,  MA  02139 

1  Niedlinger  Associates 
Consulting  Eng 
110  East  59th  Street 
New  York,  NY  10022 
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2  Virginia  Polytechnic 
Inst  i  Stats  Uni/ 

Dept  of  Eng  Mach 
ATTN:  Prof  D.  Frederick 
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1  Brown  University 
Div  of  Engineering 
ATTN:  Prof.  P.  Sy.nonds 
Providence,  RI  02912 

1  University  of  Dayton 
Research  Institute 
ATTN:  Dr.  S.  Bless 
Dayton,  On  45463 

1  University  of  Wasnington 
Dept  of  Aeronautics  and 
Astronautics 
ATTN:  Prof.  I.  M.  Fyfe 
Seattle,  'NA  33193 

1  Aeronautical  Research  Associates 
of  Princeton 
ATTN:  Dr.  J.  A.  Leech 
P.0.  Box  2229 
50  Nasnington  Road 
Princeton,  NJ  03540 

1  New  Mexico  Institute  of  Mining  £ 
Technology 
TERA  Sroup 
Socorro,  NM  37331 

1  Lockheed  Missiles  £  Space  Co. 
ATTN:  J.  J.  Murphy,  Dept.  31-11 
Bldg.  154 
P.0.  Box  504 
Sunnyvale,  CA  94035 
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1  The  Mitre  Corporation 
ATTN:  Library 

Stop  B^I  50 
P.O.  Box  203 
Bedford,  MA  01730 

1  Pacific  Sierra  Research  Corp 
ATTN:  Dr.  Harold  Brode 
1455  Cloverfield  Boulevard 
Santa  Monica,  CA  90404 

1  Physics  International  Corp 
ATTN:  Technical  Library 
2700  Merced  Street 
San  Leandro,  CA  94577 

4  Carpenter  Research  Corporation 
ATTN:  Jerry  Carpenter 
J.  3.  Lewis 
Technical  Library 
Allan  .Xuhl 
P.O.  Box  9595 
Marina  del  Rey,  CA  90291 

1  Science  Applications,  Inc. 
ATTN:  Technical  Library 
1250  Prospect  Plaza 
La  Jolla.  CA  92037 

2  Science  Applications,  Inc. 
ATTN:  Burton  3.  Chanbers 

John  Cockayne 
P.O.  Box  1301 
1710  Goodridge  Drive 
McLean,  VA  22102 

1  Science  Systems  and  Software 
ATTN:  Technical  Library 
P.O.  Box  1520 
La  Jolla,  CA  92037 


2  McDonnell  Douglas  Astronautics 
Company 

ATTN:  Rooert  'A.  Halprin 
Dr.  P.  Lewis 
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•  Holifield  National  Laboratory 
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1  Washington  State  University 
Pnysics  Department 
ATTN:  3.  R.  Fowles 
Pullman,  WA  99153 
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ATTN:  H.  N.  Aoramson 
A.  B.  Uenzel 
P.O.  Drawer  23510 
San  Antonio,  TX  73234 


Aberdeen  Proving  Sround 

Dir,  U3A.1SAA 
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ATTN:  AMSTE-TO-F 
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ATTN:  SMCCR-RSP-A 
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